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ABSTRACT 

Context. Close-in exoplanets interact with their host stars gravitationally as well as via their magnetized plasma outflows. The rich 
dynamics that arises may result in distinct observable features. 

Aims. Our objective is to study and classify the morphology of the different types of interaction that can take place between a giant 
close-in planet (a Hot Jupiter) and its host star, based on the physical parameters that characterize the system. 

Methods. We perform 3D magnetohydrodynamic numerical simulations to model the star-planet interaction, incorporating a star, a 
Hot Jupiter, and realistic stellar and planetary outflows. We explore a wide range of parameters and analyze the flow structures and 
magnetic topologies that develop. 

Results. Our study suggests the classification of star-planet interactions into four general types, based on the relative magnitudes 
of three characteristic length scales that quantify the effects of the planetary magnetic field, the planetary outflow, and the stellar 
gravitational field in the interaction region. We describe the dynamics of these interactions and the flow structures that they give rise 
to, which include bow shocks, cometary-type tails, and inspiraling accretion streams. We point out the distinguishing features of each 
of the classified cases and discuss some of their observationally relevant properties. 

Conclusions. The magnetized interactions of star-planet systems can be categorized, and their general morphologies predicted, based 
on a set of basic stellar, planetary, and orbital parameters. 

Key words. Planet-star interactions - Stars: winds, outflows - Magnetohydrodynamics (MHD) - Methods: numerical 


1. Introduction 

The exoplanet class “Hot Jupiters” consists of massive gaseous 
planets (of mass -1-10 Mj, where the subscript “J” denotes the 
planet Jupiter) that orbit their host stars on close-in trajectories 
(semi-major axis of <0.1 AU, corresponding to an orbital period 
of a few days). Among other exoplanets, the size and location of 
Hot Jupiters favors both their discovery and the collection of de¬ 
tailed data. The two leading and complementary detection tech¬ 
niques, radial velocity and transit observations, have resulted in 
a wealth of information regarding their orbital charac teri stics , 
such as obliquity and eccentricity (e.g jFabrvckv & Winnl2009l) . 
as well as their physical properties, such as surface tem pera- 
ture and atmospheric composition (e.g., [Bean et aPl2Q13h . The 
study of the phenomenology of Hot Jupiters also provides valu¬ 
able input to research on the physical conditions and processes 
in smaller, h arde r to probe, Neptune- and Earth-like bodies (see, 
e.g JSeagertl201ll for a recent reference on exoplanets). 

Hot Jupiter atmospheres are believed to be heated by pho- 
toionizing extreme-u ltraviolet (EU V) stellar irradiat i on (e . g., 
iBurrows et al.l 2000i but see also lOwen & Jacksonl 120121 as 
well as iBuzasi 120131 and IT anzall2013[ who proposed, respec¬ 
tively, X-rays for very young stars and magnetic reconnection 
for very close-in planets). The energy input could be large 
enough to induce an evaporative mass loss at a rate of - 10 10 - 
10 12 g s -1 . Such outflows have been inferred in a couple of sys¬ 
tems (HD 189733b and HD 209458b ) in which detailed obser- 
vations have been carried out (e.g., Ividal-Madiar et al.l 12004 


Ehrenreich et al. 20081 Lecavelier des Etangs et al. 

201(| 

2012t 

Bourrier et al. 2(H3l although see Ben-Jaffel f2007 

2008 

for a 


different interpretation of the data). The theory of evaporative 
planetary outflows has been developed in a number of papers, 
taking into account the thermal and ionization structure of the 
irradiated planetary atmosphere as well as tidal and magnetic ef- 


fects (e.g.. Lecavelier des Etangs et al. 

2004;[Yelle 2004 

2006; 

Tian et al. 2005t iGarcia Munoz 2003" 

Murrav-Clav et al. 

2009; 

Adams |20 ill [Trammell et al. 2011. 20141: Koskinen et al. 

2013; 

Owen & Adams 20141 and references therein). 


The host stars of close-in planets typically drive magne¬ 
tized stellar winds that are accelerated to speeds of a few hun¬ 
dred kms -1 . The properties of these winds change as the star 
evolves, f rom strong and often collimated outflo ws in young 
stars (e.g.. [Matt & Pudritzll2008l:[Matt et al.l 2 012b ) to weak and 
more isotropic winds during the main -sequence phase (e.g., 
[Matt et al1l2012al:lCohen & Prakell20l4 and references therein). 
The magnetized plasma that is driven out in these winds con¬ 
stitutes the medium into which the evaporative outflows from 
close-in expoplanets expand. 

A close-in planet and its host stars could thus interact 
through their respective outflows and/or magnetic fields. In ei¬ 
ther case, the interaction could have a potentially observable 
signature. For example, it has been suggested that there might 
be a detectable bow shock f ormed by the supersonic ste lla r wind 
in front of the planet (e.g.. Ivklotto et al.l 12010.1201 lallbL 120141: 
iLlama et al]l201 lll2013UBen-Jaffel & Ballesterll2013l) . or an ob¬ 
servable cometary-type tail associated with the swept-up plane- 
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tary wind that trails the planet in its orbit (e.g., I Schneider et afl 
119981: iMura etaD 120111: iKulow et aPl2014l) . According to these 
proposals, a bow shock could potentially contribute to the ab¬ 
sorption of stellar radiation before an ingress of a transiting 
Hot Jupiter, whereas a tail might contribute after an egress. 
It was furthermore proposed that the interaction with the stel¬ 
lar magnetosphere might lead to p hased hot spots on the stel- 


lar surface or to flaring activity (Shkolnik et al. 

20031 2004 

2005 

, 2008; Walker et al. 200| Lanza 2008, 2009 

: ICohen et al.1 

2009 

; Pillitteri et aD [2010, 2011, 120141), and that the plane- 


tary magnetic fiel d could be affected in both its topology and 


its strength (e.i 

e.. Laine et al. 2008|; (Cohen et al. 201 lbL 

20 ll 

Strugarek et al. 

20 {3). It was also suggested ( Cohen et al.l 

2010) 


that the interaction with a magnetized Hot Jupiter could lead to 
a reduction in the mass loss rate, and even more strongly in the 
angular momentum loss rate, from the host star. 


In this paper we perform 3D magnetohydrodynamic (MHD) 
simulations that aim to address the following questions: 

1. How can one classify the different types of star-planet inter¬ 
actions that arise when a Hot Jupiter orbits its host star? 

2. What are the dynamical features that develop in each case, 
and what are their expected observational signatures? 

3. What is the impact of this interaction on the planetary and 
stellar outflows? 

We classify such interactions by studying the evolution of a grid 
of numerical models. In particular, we investigate whether it is 
the planetary magnetosphere or the planetary outflow that in¬ 
tercepts the stellar plasma, whether the stellar outflow termi¬ 
nates in a bow shock, the circumstances for forming a stable 
tail along the planet’s orbit, and the conditions under which ac¬ 
cretion streams of planetary material reach the stellar surface. 
We carry out a parametric study and present the criteria that can 
be used to classify a system’s morphology and dynamics based 
on its physical parameters. We also consider some general prop¬ 
erties of the interaction regions that bear on their potential de¬ 
tectability. 

Previous work by ICohen et all (1200911201 lah has identified 
and described some of the physical processes that may take place 
in a system consisting of a magnetized star and a magnetized 
Hot Jupiter. These studies were, however, focused on the inter¬ 
pretation of two specific systems. In addition, the planetary out¬ 
flows were not explicitly included; they developed based on the 
boundary conditions applied on the surface of the planet. In con¬ 
trast, we perform a systematic study to classify the types of star- 
planet interactions, incorporating consistent planetary winds that 
are based on detailed models. 

The paper is structured as follows. Sect. [2] explains the 
methodology and presents the numerical models. Sect. [3] reports 
the results of the simulations and describes the dynamics. Sect.0 
classifies the types of interaction and considers some of their po¬ 
tential observational manifestations. Finally, Sect. [5] summarizes 
the conclusions of this work. 


2. Numerical modeling 

As discussed in [M urrav-C lav et al.l (120091 see also iTian et all 
120051: [G arcfa Munozl 120071) . evaporative Hot Jupiter outflows 
need to be modeled as a fluid (using a hydrodynamic or an MHD 
formulation) rather than as a collection of individual particles 
that undergo Jeans escape, because the plasm a typically r emains 
collisi onal beyond the flow’s critical point. iMurrav-Clav et al.l 
(2009) considered an unmagnetized outflow from a Hot Jupiter 


exposed to a UV flu x ne ar the low end of the range of values that 
we model (see Sect. m, and evaluated the collisional mean free 
path l vv = l/(nH<x P p) to proton-proton collisions in a hydrogen 
plasma of density nu and temperature T « 10 4 K (for which the 
relevant cross section is <x pp « 10 _13 [T/10 4 K] -2 cm 2 ). They in¬ 
ferred a value (^ 10 7 cm) that is ~10 4 smaller than the density 
scale height at the sonic point (which lies at a distance of ~5 
planetary radii from the center of the planet). The large value of 
this ratio indicates that the collisional approximation should be 
applicable to this problem over a wide range of circumstances. 
Similarly, our simple stellar-wind model (Sect. 12.81) is also ad¬ 
equately modeled within the hydrodynamic framework (e.g., 
[Parked I 96 0). The fact that both the planetary and the stellar 
outflows are likely magnetized sharpens these conclusions: for 
typical parameters, the proton Larmor radius = c s m v c/(eB) 
(where c s is the speed of sound, m p is the proton mass, c is the 
speed of light, e is the electron charge, and B is the magnetic 
field amplitude) is ~ 10 2 (c s /10kms _1 )(B/l G) -1 cm — much 
smaller than ^pp- Given that the region where the planetary and 
stellar outflows interact is located at a distance from the source 
of each of these winds that is comparable to the distance of its 
respective critical surface, an MHD treatment of this interaction 
should be amply justified. 


A numerical simulation that encompasses both the star and 
the planet needs to adequately resolve the sizes of both objects. 
For the case of a Hot Jupiter orbiting a solar-type host, the 
characteristic length scales of the planet and the star differ by 
roughly an order of magnitude. Since 3D numerical simulations 
are computationally demanding, the best strategy to model the 
star-planet interaction is to adopt multiple levels of refinement. 
One approach is to configure an adaptive mesh that follows the 
orbit of the planet and provides high resolution for its interior 
and environment. Another is to set up a static, multiply refined 
grid and combine it properly with a corotating frame of refer¬ 
ence, so that the location of the planet stays fixed in the highly 
resolved region. Here we adopt the latter approach, which also 
simplifies the treatment of the interior of the planet. In general, 
a static, multiply refined grid is also appropriate for inclined or¬ 
bits, as long as they are circular. 


In the following, we denote the star with and the planet 
with “o”. The simulations are performed in Cartesian coordi¬ 
nates (x, y, z), but for convenience we also use the notation of 
spherical coordinates (. R , 6 , 0). The planet is assumed to orbit 
in the x-y plane, or, equivalently, at 6 = 7t/2. The center of the 
star is placed at the origin, (x*, y*, z*) = (0, 0, 0), whereas the 
center of the planet at (x Q , y Q , Zo) = (jc 0 , 0, 0). For the initial¬ 
ization we use an extra set of coordinate systems denoted with 
primes, (x\ y', z') and (R\ 6 ', 0'), with their origins placed at 
(Vo, y Q , Zo) (i.e., centered on the planet). At each spatial point 
(v, y, z) we first calculate (v', y', z') = (x - x Q , y, z) and then 
perform a coordinate transformation to obtain (R\ O', 0'). In or¬ 
der to keep the center of the planet fixed, we adopt a corotating 
frame with Of r = O 0 rb> where Of r = Df r z is the angular velocity 
of the frame and D or b = (GM*/v 3 ) 1/2 is the Keplerian angular 
speed of the planet (with M* being the mass of the star and G the 
gravitational constant). Since M* » M 0 (where M 0 is the mass 
of the planet), we assume that M* + M Q ^ M* and hence that 
the origin of the coordinate system is effectively located at the 
center of mass. 
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2.1. The MHD equations 

The ideal MHD equations in a frame that corotates with the 
planet are: 

d -£. + v-{pv) = o, (i) 

ot 

% + 0 • V)v + -j—B x(VxB) + - VP = g + Fin, (2) 

ot Anp p 


dP 

-h v • VP + yPV • v = 0, 

ot 


dB 

— + V x (Z? x v) = 0. 
ot 


(3) 

(4) 


The variables p, P, v, and B denote the density, pressure, ve¬ 
locity, and magnetic field, respectively, and y is the polytropic 
index. The magnetic field is required to obey the solenoidal con¬ 
dition V • B = 0. Since we consider the nonrelativistic regime, p, 
P, and B have the same values in the corotating frame and in the 
lab (star’s center of mass) frame. The velocity in the lab frame 
is obtained from the expression vi a b = v + sin 6R£lf r (j). The vec¬ 
tor g = g* + go is the gravitational acceleration, and the inertial 
force that appears in the corotating frame, F m = F centr + F cor , 
has centrifugal 

Fcentr — — [^fr x (Of r x R) j = + yy) (3) 

and Coriolis 


see also iTrammell et al"l201ll). but their results r emain qualita¬ 
tively similar to those of lMurrav-Clav et all (120091) . We therefore 
consider the adopted approximation to be justified for the present 
study, especially in view of the fact that the classification crite¬ 
ria presented in Sect. 0] do not depend on the details of the wind 
models. 

We further assume, for simplicity, that the plasma consists 
of pure atomic hydrogen, and that it is fully ionized, so that the 
temperature is given by T - Pm p /2pk^ (where the factor 1/2 
represents the mean molecular weight, and &b is the Boltzmann 
constant). The strong ionization and near-isothermality assump¬ 
tions for the planetary wind are consistent with the “high UV 
flux” case (corres ponding to a young host st ar; e.g.. lRibas et al.l 
120051) presented in lMurrav-Clav et ahl (l2009h . By contrast, in the 
“low UV flux” case considered by these authors (corresponding 
to a solar-analog host), the wind temperature drops with dis¬ 
tance from the stellar surface on account of adiabatic cooling. 
To mimic this case, we also consider models with a characteris¬ 
tic temperature that is lower than 10 4 K. 

At locations where the magnetized stellar wind interacts with 
the planetary outflow or magnetosphere, field-line reconnection 
may potentially take place. This process can be realized in the 
simulations — even though resistivity is not explicitly included 
in the formulation — because numerical diffusion allows the 
magnetic field to reconnect in the presence of strong current 
sheets. We note, however, that the model setup adopted in this 
paper represents the star and the planet as aligned magnetic 
dipoles, a configuration that does not give rise to X-point struc¬ 
tures and is therefore not prone to reconnection. 


Fcor — 2(Of r x v) — 2£lf r (VyX v x y) (6) 

components. 

The temperature at the base of an EUV-induced Hot Jupite r 
outflow is found to be ~ 10 4 K (e.g.. lMurrav-Clav et al .112009 ): 
this value reflects the balance between photoioniz ation heating 
and collisionally excited radiative cooling (e.g.. ISpitzerll 1978h . 
Furthermore, in cases where the stellar UV flux is sufficiently 
high, radiative cooling continues to roughly balance the pho¬ 
toionization heating as the wind moves away from the surface, 
with the result that the tempera ture changes little along the flow 
(e.g.. iMurrav-Clav et al.ll2009h . The isothermal approximation 
was also found to apply to the inner regions of the solar wind 
(e.g JCranmer et al]l2007l) . and has t hus been adopted in p revi- 
ous simulation s of both stellar (e.g., IMatt & Pudri tz 2008) and 
planetary (e.g.. ITrammell et al1l2014l) outflows. Based on these 
considerations, we model both the stellar and the planetary out¬ 
flows in our simulations as being nearly isothermal throughout, 
characterizing th em by a p olytro pic in dex y = 1.05 as in the 
fiducial model of iMatt & Pudritzj (2 0081) . This approach enables 
us to simulate the dynamics of these winds without having to 
deal with the intricacies of a detailed calculation of the thermal 
structure of the outflow. As we show in Sect. 12.81 the adopted 
framework can be made to yield density and velocity profiles 
for the planetary wind that are compatible with the results ob¬ 
tained by e mploying the more elabo rate (albeit nonmagnetic) 
ID model of iMurrav-Clav et aTl ( 2009 ). which incorporates pho¬ 
toionization heating and radiative as well as adiabatic cooling. 
iKoskinen et al. (201 3) — who adopted a more detailed model¬ 
ing approach that included, among other factors, a consideration 
of the full spectrum of ionizing photons (rather than just a single 
rep resentative frequency) — found quantitative differences from 
the lMurrav-Clav et al.l ( 2009 ) findings (such as a slower increase 
of the ionization fraction with distance from the planet’s surface; 


2.2. Stellar and planetary parameters 

Table \T\ lists the basic model parameters that characterize the 
host star and the Jupiter-like close-in planet. The masses (M*, 
M 0 ) enter only in the gravitational field at the exterior of each 
body, whereas their radii (F*, R 0 ) define the spheres within 
which an internal boundary is applied. The temperatures cor¬ 
respond to the values at the base of the respective outflows (the 
stellar corona and the UV absorption layer in the outer planetary 
atmosphere). The value of the density of the stellar outflow, p*, is 
chosen so that the mass-loss rate be comparable to that of the so¬ 
lar wind (a few t imes 10“ 14 Mq yr -1 ), an approach adopted from 
iMatt & Pudritz ( 200 8). For the planet, the base density value is 
determined from the requirement that the simulated wind match 
high-resolution ID numerical simulations that include the rel¬ 
evant heating and cooling processes (see Sect. 12.81) . Thus, p G 
should not be regarded as the actual density at optical depth 
r ~ 1 (the physical base of the outflow) but rather as a numer¬ 
ically motivated value that gives the appropriate wind profile. 
The initial magnetic fields of the star and the planet are assumed 
to be dipolar, with equatorial surface values equal to B * and B Q , 
respectively. 

The escape speed from the stellar surface is given by v esc * = 
(2GM,/^) 1/2 , whereas the sound speed at the base of the stellar 
wind is c s * = (yF*/p*) 1/2 . A useful model parameter for de¬ 
termining the properties of the wind is T* = (1/2)(v esc */c s *) 1 2 . 
Another important nondimensional variable is the ratio of the 
thermal and magnetic pressures at the stellar surface, /?* = 
F*/(F*/87t) (the stellar plasma beta), which characterizes the dy¬ 
namical significance of the magnetic fieldQ The third relevant 
parameter for characterizing the outflow is 6* = (Q*F*/c s *) 2 , 


1 For plasma-/? values less than 1, the magnetic field dominates 

the dynamics and constrains the plasma to move along flux tubes. 
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Table 1 . Range of parameters considered in the simulations. 


Parameter 

Symbol 

Stellar value 

Planetary value 

Units 

Radius 

R*,R 0 

1% 

1-1.5 Rj 

Rj * 1CT 1 Rq 

Mass 

M*, M 0 

1 Mq 

0.5-1 Mj 

Mj ^ 10" 3 Mq 

Temperature at the base of the outflow 

T*,T 0 

10 6 

6 x 10 3 -10 4 

K 

Density at the base of the outflow 

P*,Po 

5 x 1(T 15 

7 x 10 -17 —3 x 10" 13 

gem -3 

Equatorial surface magnetic field 

B„B 0 

2 

0.1-1 

G 

Escape speed at the base of the outflow 

^esc *> r esc o 

620 

35-60 

kms -1 

Sound speed at the base of the outflow 

C s *> c so 

130 

9-13 

kms -1 

Escape velocity parameter (1/2)(v esc /c s ) 2 

4*, 4 0 

11.5 

3.8-23 

- 

Plasma [P/(B 2 /8n)] at the base of the outflow 

/5„/5o 

5 

0.002-400 

- 

Rotation period 

P„Po 

12 

1.2-3.7 

days 

Rotation speed parameter (£lR/c s ) 2 

6*, e 0 

0.1 

0.03-0.5 

- 

Lagrange point 1 

Lu, L\ q 

4.6-9.5 Rq 

2.7-6.9 R] 

Rj * 1CT 1 Rq 

Stellar UV flux 

F V v 

5 x 10- 5 x 1CP 

- 

erg cm -2 s -1 

Orbital radius 

Rorb 

- 

0.023-0.047 

AU 

Orbital period 

n>rb 

- 

1.2-3.7 

days 

Orbital speed 

Vorb 

- 

138-195 

kms -1 


which measures the effect of surface rotation. Corresponding ex¬ 
pressions for v esC o, c so ,do,/?o, and e Q are used for the planet. 

The stellar UV flux, F\jy, is a critical quantity because it 
determines the energy input and hence the strength of the plan¬ 
etary wind. Finally, Table |T| gives the rotational periods of the 
star (!P* = 27r/Q*) and the planet (P Q = 27r/Q 0 ) as well as 
the orbital characteristics, * or b = x Q , P Q vb - 27r/Q 0 rb, and 
v 0 rb = (GM*/R 0 rb) 1/2 . The host star is taken to rotate at 1% of 
the breakup angular speed (as compared to 0.4% of ( GM*/Rl ) l/2 
for the present sun), whereas the planet is assumed to be tidally 
locked (i.e., O 0 = n or b) ? which is a reasonable approximation 
for a close-in planet tha t interacts tidally with its host star (e.g., 
iMatsumura et al.l [20 1 Oh . Our model planets thus rotate at ~3- 
25% of breakup. 

2.3. Stellar and planetary winds 

We base the numerical treatme nt of the two outflows on the sim¬ 
plified approach of iMatt & Balick ( 2004 ). In brief, we initialize 
dipolar magnetospheres (force-free by definition) and keep fixed 
at each surface the physical quantities that drive the flow. The 
temporal evolution of such configurations opens u p the magne¬ 
tospheres and leads to steady-state winds (see IMatt & Pudritzl 
120081) . Although this method does not require the winds to be 
specified explicitly, we have opted to initialize a simple form of 
outflow to ensure that the starting wind profiles are in agreement 
with detailed, self-consistent models in the literature. We em¬ 
phasize that, because of the transonic nature of these flows, the 
final steady state does not depend on the initial configurations, 
only on the boundary conditions that are imposed at the surfaces 
of the two objects. However, in low-resolution simulations that 
do not include all the relevant physics, simple assumptions about 
the density and temperature at the boundary might not give cor¬ 
rect results. Therefore, we use the profiles of detailed outflow 
models as a guide in order to set up the appropriate boundary 
conditions. This procedure is particularly useful for capturing 
the correct mass loss rate of the planet, as described in more de¬ 
tail in Sect. 12.81 

We initialize the velocity field of each out flow a t time t - 0 
as an isotropic, isothermal Parker wind ([Parked 19581) . In particu- 


Conversely, when the value of this parameter exceeds 1, the fluid dom¬ 
inates the dynamics and the magnetic field is dragged by the flow. 


lar, the velocity profile for the stellar outflow, v™*(*), is obtained 
by solving numerically the equation: 

^-ln^ = -3-41n(y) + lnf. + 2p (7) 

where i//*(R) = (v^*/c s *) 2 > £*(*) = */**, and 4* is defined 
in Sect. 12.21 In this equation we take the speed of sound to be 
strictly isothermal, i.e., c s = (2k^T/m v ) l/1 (corresponding to a 
poly tropic index that is exactly 1). Since the star is assumed to be 
rotating, we make the approximation that the initial wind is ro¬ 
tating rigidly with it. This is a minor contribution to the outflow 
velocity because Q* is an order of magnitude smaller than Qf r 
(and, in any case, the flow self-consistently assumes the correct 
values of as the simulation evolves). The Cartesian compo¬ 
nents of the initial stellar outflow velocity are then given by the 
sum of the wind, rotation, and frame speeds: 


v'^(x, y, z) = Sin 0 1 

cos <b v™‘<7?) + sin 0 R(Q b + Q„)] x, 

(8) 

v'yl\x,y,z) = sin (91 

sin 4> v^'l(R) - cos 4> R(Q fl + Q*)] y, 

(9) 

v“ t (x,);,z) = cos0v“(«)z. 

(10) 

To derive the corresponding expressions for the planetary out¬ 
flow, we first calculate v^C*') from Eq. Q, using if/ 0 (R '), 
and (*') = *7*o. We then add the planet’s orbital velocity, 
v 0 rb = *orb^orbJ ? as well as the frame and planetary-rotation ve¬ 
locities (using the rigid-body assumption), to obtain: 

vfXx,y,z) = sin#' 

[cos 



+ sin (f> RQ. b - sin <p'R'Q 0 J x, 

(ii) 

v^(x,y,z) = sine' 

[sin tp'v™l(R') 



- cos 0 RDf r +R 0 rb^orb + cos 0'R'Do] y , (12) 

v'^(x,y,z) = cos e'v^ 0 (R')z, 

(13) 


where *', O', and 0' are functions of (x, y, z), as explained at 
the beginning of this section. For simplicity, we have assumed 
that the initial planetary wind orbits and spins rigidly with the 
planet, an approximation that is appropriate only close to the 
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surface. Note that, for a tidally locked planet, the orbital veloc¬ 
ity cancels out with the frame and planetary rotation velocities, 
and the planet stays fixed at the location (x G , y Q9 z Q ). For exam¬ 
ple, the y components of these three terms at a distance v along 
the line that connects the centers of the star and the planet are 
-vQf r + x o n o rb + v'Q 0 = 0. Under this assumption, given that we 
are working in the corotating frame, we could have just as well 
initialized the velocity by using only the thermal (Parker) wind, 
with no extra terms. 

To obtain the pressure profile of the stellar wind, we solve 
analytically the hydrodynamic radial momentum equation, 


dvR dP GM * 

PVR ~dR + dR ~ P R- 
which gives: 


P‘* nit (x, y, z) = P* exp 



(14) 


(15) 


The density is then found using the isothermal-plasma assump- 
tion: 


pi ni w,z) = 



(16) 


In a similar vein, we get for the planet 


PT(x,y,z) = P 0 exp 



and 


pl nit (w) 



(17) 


(18) 


We set up the stellar wind everywhere in the computa¬ 
tional box, apart from a sphere of radius 10 R 0 , centered on the 
planet. Within this volume we initialize the planetary wind. Both 
outflows are initially supersonic and super-Alfvenic (with the 
Alfven speed given by va = [B 2 / 4np] l/2 ) at the interface that 
separates themQ 

It is worth noting that the Parker wind (Eq. [7]) is nondimen- 
sional (with the radial distance expressed in units of the radius of 
the object and the velocity in units of the sound speed at the base) 
and that its solution depends only on the parameter A, which can 
be written in the form 


A 


1 / 2GM/R \ = /JM/M/M 

2 \2k B T/mJ R T J ' 


(19) 


Interestingly, for a solar analog ( R * = 1 Rq, M* = 1 M©, 
T* = 10 6 K) and a Hot Jupiter (. R 0 = IRj - 10 ~ l Rq, M 0 = 
1 Mj ^ 10 _3 M©, T 0 = 10 4 K) we obtain A* « A 0 . As a re¬ 
sult, the two outflows reach their corresponding sonic speeds at 
the same scaled distance. For the adopted fiducial parameters, 
A = 11.5, which implies that the stellar and planetary outflows 


2 Recall that in MHD there are three types of waves (listed in order 
of increasing propagation speed): slow-magnetosonic (SMS), Alfven, 
and fast-magnetosonic (FMS). Along field lines, an FMS wave propa¬ 
gates at max(v A , c s ) and an SMS one at min(v A , c s ). Across field lines, 
the SMS and Alfven wave propagation speeds are zero, and an FMS 
wave propagates at (v^ + c 2 ) 1/2 . The FMS critical surface is the separa- 
trix beyond which no perturbation can propagate oppositely to the flow 
(analogous to the sonic critical surface in hydrodynamics); we hence¬ 
forth refer to this surface as the “critical surface”. We also follow a com¬ 
mon practice in the literature and plot the Alfven and sonic separatrices 
rather than the SMS, Alfven, and FMS ones. 


become supersonic at a radius that is slightly smaller than 6 R* 
and 6R 0 , respectively. 

In a real system, the stellar flux that heats the planetary 
surface varies with longitude (in particular, there is no di¬ 
rect irradiation of the planet’s night side) as well as latitude 
(in particular, the planetary poles also receive no direct heat¬ 
ing). This implies that even an unmagnetized outflow would be 
anisotropic, although the quantitative effect of the uneven irradi- 
ation_on the evaporative mass outflow rate need not be large (e.g., 
iMurrav-Clav et al.ll2009l) . In this work we simplify the treatment 
by assuming that the base temperature and density are uniform 
over the entire planetary surface. 

The initial magnetic field configuration in the computational 
domain is taken to be a sum of two magnetic dipoles, + 
which are given by 

B R 2 

Bf\x,y, z) = [3xzx + 3yzy + (3r - R 2 ) z] (20) 

and 

Bf\x, y,z)= 20- [3.vY X + 3vY y + (3z' 2 - R' 2 ) z] (21) 

for the star and planet, respectively. Both dipoles are oriented 
along the z axis and are aligned with each other. We also write 
explicitly the components of the total gravitational field g out¬ 
side the two bodies: 

GM * GM * 

g*(x,y,z) = -fir R = -p-(xx + yy + zz) , (22) 

go(x, y, z) = (x x + y' y + z' z) . (23) 


2.4. Wind boundary conditions 

In order to continuously accelerate the outflows and guarantee 
a constant supply of mass, we keep fixed the density, pressure, 
and velocity at their bases. Specifically, during the temporal evo¬ 
lution, we impose the initial values of these quantities at the re¬ 
gions defined by R* < R < 1.5 R* and R Q < R' < 1.5 R Q for 
the stellar and planetary winds, respectively. However, the mag¬ 
netic field is free to evolve within these shells, readjusting self- 
consistently from its initial dipolar topology. This is_a simplified 
treatment as compared with the formulation of iMatt & Pudritzl 
( 2008 ). who modeled 2.5D axisymmetric stellar winds using a 
thinner, four-layer shell above the surface of the star@ Within 
that region, they progressively relaxed the time-dependency con¬ 
straint on the physical quantities. Here we employ a less de¬ 
tailed setup on account of the lower resolution imposed by our 
3D modeling. Nevertheless, the steady-state winds that we ob¬ 
tain compare well with the more refined 2.5D simulations (see 
Sect. ESJ). 

The above implementa tion of stell a r wind s d iffers from 
the approach adopt ed by ICohen et all (1201 lallbL 120141) and 
lYidotto et al.l (120141) (and references therein), which includes 
nonaxisymmetric and/or temporal variabilities. We neglect such 
effects in this work in order to investigate the star-planet inter¬ 
action in a general manner that does not depend on temporary 
local fluctuations. 

Our wind boundary conditions do not incorporate the ef¬ 
fect of tidal forces on the planetary outflow. As discussed by 

3 A 2.5D MHD simulation refers to the evolution of all three compo¬ 

nents of the velocity and magnetic field in a 2D computational domain. 
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[Trammell et al.ll ( 20111 these forces could suppress the wind 
from the polar regions of a tidally locked planet with a dipolar 
field geometry. In particular, these authors showed that the wind 
would not undergo a sonic transition in this case if e Q > 4/T 2 
(where the parameters e and A are defined in Sect. 12.21) . Our ne¬ 
glect of this suppression in the adopted boundary conditions is 
not a serious omission in view of the fact that our simulations 
incorporate the effect of these forces for 7?' > 1.5 R 0 and that 
the sonic surface typically lies at a distance of a few R 0 from 
the planet’s centerQ Furthermore, for the star-planet configura¬ 
tions considered in this paper, which are characterized by paral¬ 
lel orbital and spin axes, the bulk of the interaction between the 
stellar and the planetary plasmas occurs in the equatorial plane, 
with the details of the planetary outflow from the polar regions 
having little effect on the results. 

2.5. Stellar and planetary interiors 

The interiors of both the star and the planet are treated as an in¬ 
ternal boundary. In particular, the values of all physical quanti¬ 
ties are kept fixed by having them overwritten at each time step. 
To avoid spurious effects at the surfaces of the objects, or ex¬ 
treme dynamics that could affect the global time step of the sim¬ 
ulation, we prescribe an approximate equilibrium within their 
volumes. 

Since the physical conditions in the interior of the star or 
the planet are not important for the simulation, the simplest ap¬ 
proach is to consider them as uniform-density spheres. The inte¬ 
rior gravitational accelerations are then given by 


g*(x,y, z ) = -nGp t R = -nGp*(xx + y y + zz) 

(24) 

and 


4 

go(x, y, z) = -nGp 0 (x’x + y'y + z z ), 

(25) 

from which the pressures can be inferred using the hydrostatic 
equilibrium condition: 

Pf\x,y,z) = P, + \nGp\ [r] - R 2 ) , 

(26) 

Pf\x,y,z) = Po + 2 -nGpl (r 2 c - R' 2 ) . 

(27) 


We keep the dipolar expression for the magnetic field, but 
we smooth it out at the center where it becomes singular. 
Specifically, we assume that the central regions of the star and 
the planet, 7? < 7?*/2 and R' < RJ 2, are uniformly magnetized 
spheres, with magnetic fields Bf lt = 167?*z and Bf 11 = 16 B 0 z 
that connect smoothly with the dipolar configurations at larger 
radii. Finally, we approximate the star and the planet as rotating 
rigid bodies and specify their angular velocities by O* = Q*z 
and Q 0 = Q 0 z, respectively. 

2.6. The models 

Table [2] lists the models that we study numerically. We aim to 
explore the observationally relevant regions of parameter space 
and attempt to deduce, based on the behaviors exhibited by the 

4 We note in this connection that the only two models in Table [2|that 
develop supercritical outflows, FvRb and FvRB, do not satisfy the in¬ 
equality £ 0 > 4 /A 2 q . This is consistent with the fact that the sonic surface 
plotted in Fig.[9]is represented by a closed curve. 


simulated systems, a general — and largely model independent 
— classification scheme for magnetized star-planet interactions. 
We usually adopt two representative values for each of the phys¬ 
ical variables that need to be specified in the model. Thus, we 
consider two combinations of planetary parameters, the first hav¬ 
ing Jupiter’s mass and radius and the other corresponding to a 
less massive (0.5 Mj) and larger-radius (1.5 7?j) planet. These 
choices are intended to sample the typical range of values for 
the escape speed v esC o °c ( M 0 /R 0 ) l/ 2 from a hot Jupiter, which is 
relevant to the mass evapo ration rate M ev in the low -TAv limit 
(where M ev 00 l/Ve SCO ; e.g jMurrav-Clav et al.ll2009h . 

We similarly consider two magnitudes for the UV flux, 
one low (Fuv = 5 x 10 2 erg cm -2 s -1 ) and one high (Tqjv = 
5 x 10 5 erg cm -2 s -1 ), adoptin g numerical values s i milar to those 
given for these two limits in iMurrav-Clav et all (120091) . These 
two values correspond to different evaporation regimes. In the 
high-flux case, the excess photoionization energy is largely lost 
by radiative cooling (radiation/recombination-limited evapora¬ 
tion), whereas in the low-flux case it is mostly absorbed and 
drives the outflow (energy-limited evaporation). The outflow 
from a strongly irradiated planet is denser, hotte r, and faster 
than that from a weakly irradiated one (e.g.. IMurrav-Clav et al] 
l2Q09h . Since in this paper we do not explicitly model the radia¬ 
tive processes in the planetary atmosphere, we specify the wind 
properties through the values of the base temperature and den¬ 
sity, which are adjusted to match the results of more detailed 
calculations (see Sect. 12.8k 

We also vary the distance of the Hot Jupiter from the host 
star. Smaller orbital radii imply a stronger tidal force as well 
as a locally weaker, and possibly not fully developed, stellar 
wind0 The value of the orbital radius also determines whether 
the planet lies within the critical surface of the stellar wind or 
whether it is impacted by a super-critical flow. 

The two values that we adopt for the magnetic field ampli¬ 
tude at the planet’s surface differ by an order of magnitude and 
are intended to represent a “strong” and a “weak” field. Stronger 
fields correspond to higher magnetic pressure and tension and 
exhibit a greater “rigidity.” In this case the field resists being 
opened up by the stellar outflow over a larger region near the 
equator, and as a result the planetary magnetosphere presents a 
larger obstacle in its interaction with the stellar wind. 

Although we vary the stellar UV flux, our simulations em¬ 
ploy only a single set of stellar wind parameters. Our coverage 
of the relevant parameter phase space is, however, sufficiently 
broad to yield general results, which in principle apply also to 
systems with either a weaker or a stronger stellar wind. 

2.7. Numerical setup 

The simulations are performed with PLUTO (version 4.0.1), 
a cod e for computational astrophysics dMignone et al.l I2007L 
l2012l) Fl The computational domain consists of a cube, with 
x,y,z e [-15, 15]7?*, which is resolved by a static, multiply 
refined grid of 224 x 192 x 192. We have also carried out one 
simulation (model FvRb) in higher resolution, 424 x 320 x 320, 
in order to validate our results. The region around the star, 
x,y,z £ [-5, 5]7?*, is resolved with Ax, Ay, Az ^ 0.167?* (64 3 
cells) in most cases, and with Ax, Ay, Az ^ 0.087?* (128 3 cells) 


5 We do not take into account changes to the incident UV flux that re¬ 
sult from variations in the orbital radius since their magnitudes remain 
small in comparison with the difference between our chosen represen¬ 
tative values. 

6 PLUTO is freely available at http: //plutocode. ph. uni to. it 
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Table 2. Parameters of the numerical models. From left to right: orbital radius (7? 0 rbX planetary radius ( R Q ) and mass (M 0 ), surface 
escape speed (v esC o), distance to the L\ point (in units of R 0 ), density and (p Q ) temperature (T 0 ) at the base of the outflow, equatorial 
surface magnetic field (B 0 ), and irradiating flux (high or low; see text for the numerical values). The labels of the models are divided 
into four parts: the first denotes the stellar UV flux (high [F] or low [f]), the next describes the magnitude of the escape speed (large 
[V] or small [v]), then its distance from the host star (far [R] or near [r]), and finally the strength of its magnetic field (strong [B] or 
weak [b]). 


Model 

florb (AU) 

Ro (R,) 

M 0 (Mj) 

v esC o (km s -1 ) 

L\ 0 (R 0 ) 

P„ (g cm -3 ) 

To (K) 

Bo ( G ) 

Irradiation 

FvRB 

0.047 

1.5 

0.5 

35 

3.7 

7 x 10- lb 

10 4 

1 

High UV flux 

FvRb 

0.047 

1.5 

0.5 

35 

3.7 

7 x lO" 16 

10 4 

0.1 

High UV flux 

FvrB 

0.023 

1.5 

0.5 

35 

1.8 

7 x If)" 16 

10 4 

1 

High UV flux 

fvRB 

0.047 

1.5 

0.5 

35 

3.7 

7 x 10“' 7 

6x 10 3 

1 

Low UV flux 

fvRb 

0.047 

1.5 

0.5 

35 

3.7 

7 x lO" 17 

6x 10 3 

0.1 

Low UV flux 

fvrB 

0.023 

1.5 

0.5 

35 

1.8 

7 x lO" 17 

6x 10 3 

1 

Low UV flux 

FVRB 

0.047 

1 

1 

60 

6.9 

10" 13 

10 4 

1 

High UV flux 

FVRb 

0.047 

1 

1 

60 

6.9 

io - 13 

10 4 

0.1 

High UV flux 

FVrB 

0.023 

1 

1 

60 

3.4 

10“ 13 

10 4 

1 

High UV flux 

fVRB 

0.047 

1 

1 

60 

6.9 

3 x 10" 13 

6 x 10 1 

1 

Low UV flux 

fVRb 

0.047 

1 

1 

60 

6.9 

3 x 10“ 13 

6x 10 3 

0.1 

Low UV flux 

fVrB 

0.023 

1 

1 

60 

3.4 

3 x IO" 13 

6x 10 3 

1 

Low UV flux 


in the high-resolution simulation. However, for the Hot Jupiter 
and its surroundings, x',y',z' € [-5, 5 ]R Q , we use a resolu¬ 
tion that is higher by one order of magnitude, i.e., Ax, Ay, Az - 
0.016 R* (64 3 cells) and Ax, Ay, A z - 0.008 R* (128 3 cells) for 
the standard and highly resolved cases, respectively. In the re¬ 
gion between the star and the planet we apply the resolution 
Ax, Ay, Az - 0.167?*. 

The total time of the simulation is typically ~4 days, only a 
fraction of which is required for the system to attain a quasi¬ 
steady state. For comparison, a stellar wind propagating at 
300kms _1 travels ~10 times the distance from the star to the 
outer boundary over this time interval. We adopt highly accurate 
numerical schemes to compensate for the resolution limits im¬ 
posed by the three-dimensional character of the simulations. In 
particular, we use the HLLD Riemmann solver and apply third- 
order accuracy in space (piecewise parabolic interpolation) and 
time (3 rd order Runge-Kutta). The condition V • B = 0 is im¬ 
plemented using the divergence-cleaning method, an approach 
based on the generalized Lagrange multiplier (GLM) formula¬ 
tion (see lMignone et al.ll20 Inl and references therein). 

2.8. Verification and calibration of 3D wind models 

The acceleration and final steady-state properties of a simulated 
wind depend on the included physics as well as on the resolu¬ 
tion of the grid. 3D simulations are computationally expensive 
and therefore cannot incorporate the same level of detail that 
is possible in ID and 2D models. In this work we have imple¬ 
mented a simplified procedure for treating the stellar and plane¬ 
tary outflows, and we now check on the validity of the adopted 
approximations. We first compare our 3D stellar-wind model 
with a more detailed (and higher-resolution) 2.5D model and 
verify its consistency. We then demonstrate our method of en¬ 
suring the consistency of 3D planetary wind models with ID 
hydrodynamic calculations that incorporate the relevant physics 
of evaporative outflows. 

Figure [I] compares the 3D stellar wind model adopted in this 
work (left panel), with a 2.5D simulation with the same param¬ 
eters (right panel), but with a 10 time_s_higher reso lution along 
each direction (modeled as in iMatt & Pudritzl [20081) . The den¬ 
sity, velocity, and magnetic field have very similar profiles, de¬ 
spite the fact that every cell on the left would be resolved by a 


3D stellar wind 2D stellar wind 



02468 10 02468 10 

x (R.) x (R.) 


Fig. 1. Logarithmic density contours, field lines (solid lines), 
and poloidal Alfven surface (dashed line) in the final steady state 
of the 3D stellar wind model employed in this paper (left panel) 
and of a 2.5D high-resolution simulation (right panel). 




Fig. 2. Logarithmic density (left) and radial velocity (right) in 
the final steady state of the 3D stellar wind model as a function 
of the distance R from the center of the star along the equatorial 
(v) and axial (z) directions (solid and dashed lines, respectively). 


thousand cells on the right if it were a 3D setup. Any discrep¬ 
ancies can be neglected since the outflow serves its purpose of 
providing a generic stellar wind. 

Figure [2] shows the profiles of the density (left panel) and 
radial velocity (right panel) for the steady state reached by our 
stellar wind model. The radial velocity is zero along the equator 
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Fig. 3. Logarithmic density (left) and radial velocity (right) as 
a function of distance from the planet’s center for the isotropic 
outflow model used to initialize the 3D wind simulation for a 
“light” Hot Jupiter (M 0 = 0.5 Mj, R Q = 1.5 Rj). The symbols 
“+” and “x” denote, respectively, the high- and low-flux cases 
(modeled with T 0 = 10 4 K and T 0 = 6 x 10 3 K, respectively). 
The solid and dashed lines are from the corresponding ID sim¬ 
ulations that take into account the relevant heating and cooling 
processes (see text for details). 




Fig. 4. Same as Fig. [3] but for a Jupiter-like planet (M 0 = Mj, 
Pc = Rj). 


up to ~4P*, corresponding to the extent of the stellar dead zone 
(see Fig.Q]). More generally, the finding that the wind speed re¬ 
mains lower along v than along z for any given value of R can 
be attributed to the fact that a smaller fraction of the field lines 
open up near the equator than near the poles. This behavior is 
also exhibited by more d etai led m odels of winds from solar-type 
stars (e.g. JCohen & Prakel20l4) . In the case of faster rotators, in 
which magnetocentrifugal effects play a more prominent role in 
the acceleration of the wind, the outflow speed along th e equator 
can become higher than along the spin axis (e.g jMatt & Pudritzl 
120081) . This regime is pertinent to our models of tidally locked 
planets. 

To ensure that the planetary winds are properly initialized, 
we specify the boundary conditions so that the resulting pro¬ 
files match the results of a more elaborate (and highly re¬ 
solved) ID model that incorporates relevant additional physics. 
Specifically, we simulate the outflow along the line joining the 
centers of the planet and the star, incorporating explicit energy 
and ionization balance e quations following the formulation of 
iMurrav-Clav et al.l ( 2009 ). The energy equation includes advec- 
tive and PdV -work terms as well as source terms in the form 
of photoionization heating (assuming a single ionizing photon 
energy of 20 eV and perfect efficiency) and Ly a cooling (un¬ 
der the assumption of collisional excitation by free electrons). 
This equation is integrated using a modified version of the 


Simplifi ed Non- Equili brium Cooling (SNEq) module of PLUTO 
dTesileanu et al.1112008 ). The hydrogen ionization balance equa¬ 
tion accounts for photoionization, Case B radiative recombina¬ 
tion, and ion advection. Using the results of these calculations, 
we fine-tuned the boundary conditions of the simplified 3D mod¬ 
els until the latter were found to effectively capture the general 
features of the ID models for our representative planets in both 
the high- and low-flux limits. The steady-state planetary wind 
profiles obtained in this way are shown in Figs. [3] and 0] together 
with the corresponding ID results. 

As can be seen from Figs. [3] and 01 the ID outflow models 
predict a very steep density drop near the surface of the planet. 
This implies that, if one were to choose the value of the surface 
density p D for the 3D simulation simply on the basis of the loca¬ 
tion of the nominal base of the flow (where r ~ 1), this would 
likely not lead to a good match with the ID density profile. On 
the other hand, if one were to choose even a slightly different 
reference radius for fixing the surface density, this might result 
in a significant under- or over-estimate of the mass outflow rate. 
These considerations strongly suggest that the correct procedure 
for modeling the underlying planetary outflow is to choose the 
boundary values of density and temperature on the surface of the 
planet so that they lead to a good match with the reference pro¬ 
files. The small discrepancy seen in the velocity profile of the 
low-flux case in Fig. 0] can be ignored because, as we show be¬ 
low, the interaction with the stellar plasma chokes the outflow 
from a Jupiter-like planet for this value of the flux. 

3. Numerical results 

3.1. Unmagnetized star-planet interactions 

Before turning to the MHD simulations, we make a few quali¬ 
tative remarks on the expected behavior in the absence of mag¬ 
netic fields. Assuming that there is no planetary outflow, the stel¬ 
lar wind will impact the planet directly. For a supersonic stellar 
flow, the height above the planet’s surface where the plasma will 
get shocked is determined by the location where the atmospheric 
thermal pressure, Pth o, is equal to the local value of the wind ram 
pressure, p w *v^*- For a subsonic flow, there will be no shock, 
but rather a smooth transition of the physical quantities between 
the two bodies: the density and pressure will have high values 
at the surface of each object and will be lower in between. The 
converse will hold for the velocity: zero values at the surfaces, 
and a smooth acceleration — and subsequent deceleration — on 
moving from the host to the Hot Jupiter. 

If a planetary outflow is present, its strength will determine 
whether it will be suppressed by the stellar wind (weak out¬ 
flow) or else become supersonic and then collide with the stel¬ 
lar plasma well away from the planet’s surface (strong outflow). 
When both flows are supersonic, they interact at the location 
where the ram pressure of the planetary wind, p w *v^*, is equal 
to that of the stellar gas (as measured in the corotating frame). 

3.2. General behavior 

We start our presentation of the MHD simulations by describing 
the general structure of the interacting outflows; in the ensuing 
subsections we present the results for each of the simulated mod¬ 
els in greater detail. We use model FVRB for this illustration, and 
show its initial configuration in Fig.[5j the initial structure of the 
other models is similar except for the planetary wind profiles, 
which depend on the choice of surface parameters (see Figs. 0] 
and 0]). The star is located at the center of the left panel and of 
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orbital plane face on 



-6 -4 -2 0 2 4 6 




4.2 4.4 4.6 4.8 


orbital plane face on 


orbital plane edge^on 


orbital plane edge or 


0.2 


0.0 


4.2 4.4 4.6 4.8 5.0 


Fig. 5. Initial conditions in the lab frame for model FVRB. The 
axes are in units of 0.01 AU. Logarithmic density contours (in 
units of gem -3 ) are shown in the x-y plane (z = 0; left panel) 
and in the x-z plane (y = 0; right panels). The top right panel is 
focused on the star and the bottom right on the planet. Solid lines 
represent the magnetic field, arrows show the velocity field (with 
the longest one corresponding to ~290km s _1 ), dashed lines the 
poloidal Alfven surface (y 2 + v\ — (B 2 X + B 2 )/4np ), and dotted 
lines the poloidal sonic surface (v 2 + v 2 = yP/p). The cross in 
the bottom right panel denotes the location of the L\ point. 


the top right panel, which respectively provide face-on and edge- 
on views of the orbital plane. The stellar wind is depicted with 
arrows, and its initially dipolar magnetic field with solid lines. 
Next to the star, at x ~ 4.5 (in units of 0.01 AU), is the spheri¬ 
cal volume within which the Hot Jupiter and its outflow are lo¬ 
cated (left and top right panels). A close-up of the region near 
the planet is shown in the bottom right panel. Note that both the 
stellar and the plentary winds are super-Alfvenic and supersonic 
at the initial interface between them (a sphere of radius 10F o , 
centered on the planet). 

At the beginning of the simulation, the plasma escaping from 
the surface of each object forces the magnetosphere to open. In 
general, depending on the strength of the magnetic field and the 
physical conditions at the base of the wind, this might happen 
over the entire surface (for a weaker field and/or a stronger out¬ 
flow) or just at the polar regions (for a stronger field and/or 
a weaker outflow). For the adopted value of the stellar sur¬ 
face magnetic field (. B * = 2 G), the star forms a dead zone for 
7r/4 < 0 < 37r/4. In this region, the hot plasma cannot overcome 
the opposing forces of stellar gravity and magnetic tension, both 
of which are almost perpendicular to the flow. The behavior of 
the planetary magnetosphere will be discussed on a case by case 
basis: its structure depends both on B 0 and the planetary wind 
properties and on its environment. 

3.3. Models exhibiting a planetary tail 

Figure [6] displays the final steady state of model fvRb. This case 
corresponds to a low UV flux, and the planetary outflow is weak 
and cannot overcome the ram pressure of the stellar wind (pro¬ 
jected along the orbit). The stellar wind is stopped by the mag¬ 
netic pressure of the planet’s closed field lines (which exceeds 
the thermal pressure above the planet’s surface), terminating in a 
bow shock that stands off a few planetary radii from the planet’s 
surface. The magnetosphere is dragged backward by the stellar 
wind, forming a cometary-type tail. This process is akin to the 
formation of the geomagnetic tail in the interaction of the so¬ 
lar wind with Earth’s magnetosphere (e.g jNes si ll 9651) . The tail 
does not follow the trajectory of the planet since it is continu- 


Fig. 6. Logarthimic density (color contours in units of g cm -3 ), 
velocity (arrows), magnetic field (solid lines), critical surfaces 
(dotted line: sonic, dashed line: Alfven), and the L\ point (cross) 
for the final steady state of model fvRb. Distances are measured 
in units of 0.01 AU. The panels are arranged as in Fig. 0 with 
the left and top right panels centered on the star, and the bottom 
right one on the planet. 


fvRB FVRb 



-6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6 


Fig. 7. Same as the left panel of Fig. 0 but for models fvRB 
(left) and FVRb (right). 


ously pushed away radially by the stellar wind, resulting in a 
tilted structure. 

Figure |7] illustrates how the structure of the tail is modified 
when the model parameters are changed. The left panel shows 
the results for model fvRB, which has the same parameters as 
model fvRb (Fig. [6j except that the planetary magnetic field is 
a factor of 10 stronger, corresponding to an increase by a fac¬ 
tor of 100 in the magnetic pressure at a given distance from the 
planet. As expected, this results in the bow shock being located 
farther away from the planetary surface, which, in turn, gives 
rise to a wider tail. The right panel presents the final state for 
model FVRb, which corresponds to a higher UV flux (which pro¬ 
motes an outflow) but also a larger escape speed (which hinders 
a wind). Overall, its behavior is very similar to that of model 
fvRb, although the tail in this case is noticeably denser. This can 
be understood from the fact that, in the high-flux limit, the base 
density scales roughy as F^y dMurrav-Clav et al1l2009h . 

The snapshot shown in the right panel of Fig. [7] highlights 
a few interesting dynamical features of cometary-type tails. The 
plasma that trails the planet has a velocity comparable to the 
orbital speed, whereas the component of the stellar wind along 
the orbit is negligible at that location. This results in the devel¬ 
opment of strong shear, which may trigger a Kelvin-Helmholtz 
instability. Furthermore, the stellar plasma that pushes the tail 
outward is of lower density than the trailing stream, which can 
induce a Rayleigh-Taylor instability. In this particular simula¬ 
tion, their effect on the tail structure remains comparatively mild, 
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Fig. 8. Same as Fig. [6] but for model fVRb. 



Fig. 10. 3D representation of the density structure and of the 
magnetic field lines (blue: stellar; red: planetary) for model 
FvRb. The logarithmic density scale is chosen for visualization 
purposes and does not correspond to that of Fig. 0 


likely on account of the relatively high density of the trailing 
stream. We have, however, found that the induced distortions can 
be more persistent in other cases. 

When the UV flux remains low but the escape speed is high, 
the planetary outflow weakens and the density in the tail be¬ 
comes measurably lower than in the models considered so far. 
This result follows directly from the expression for the mass 
evaporation rate in the low-flux (energ y-limited) limit, which im - 
plies that M ev ^ ^uv/v^c ( see Eq. 9 in lMurrav-Clav et al.12009b . 
In this case the tail becomes thin and barely noticeable. This is 
demonstrated in Fig. [8] for model fVRb — the outflow is even 
weaker when the magnetic field amplitude is increased (model 
fVRB, not plotted). 

3.4. Models exhibiting both a tail and an accretion stream 

Figures 0 and [TO] display the final quasi-steady state of model 
FvRb. This model differs from the reference case of Sect. f33l 
(model fvRb, shown in Fig. [6]) in having a high (rather than low) 
incident UV flux. This results in the formation of a strong out¬ 
flow, which opens up the magnetosphere and becomes super¬ 
sonic (bottom right panel) 0 

The planetary outflow propagates unperturbed over a few 
planetary radii, giving rise to the eye-shaped region around the 
planet at 4.0 < v < 5.5 (left and top right panels of Fig. 0. 

7 The comparatively weak magnetic field in this case implies that 
the Alfven critical surface is located closer to the planet than the sonic 
surface. 
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Fig. 11. Same as Fig. [6] but for model FVRB. 
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Beyond that region, the planetary outflow collides with the stel¬ 
lar wind and forms a shock. This happens on both the day and 
night sides of the Hot Jupiter on account of the high orbital speed 
(which shifts the apex of the bow shock away from the substellar 
point) and the axisymmetry assumption for the planetary wind. 

A noteworthy feature of the interaction in this case is the ac¬ 
cretion of part of the shocked planetary outflow onto the host 
star. The infalling material spirals for a quarter of an orbit and 
then impacts the stellar surface. The fact that this gas falls in, 
rather than form a torus-like structure, can be attributed to the 
action of the shock and of the stellar wind, which slow down and 
disrupt the flow, and to the increase in the gravitational pull of 
the star as the gas spirals in. The rest of the outflowing planetary 
gas is, however, pushed back and forms a tail, as in the cases con¬ 
sidered in the preceding subsection. Both the inspiraling and the 
trailing streams are affected by the velocity shear with the low- 
density stellar gas and by the outward acceleration that this com¬ 
ponent induces, which trigger Kelvin-Helmholtz and Rayleigh- 
Taylor instabilities. These effects lead to the destruction of the 
tail in this case, with its fragments being pushed outward by the 
ram pressure of the stellar wind (left panel of Fig. 0. A very 
similar outcome is produced for model FvRB (not plotted), in 
which, however, the stronger planetary magnetic field results in 
the Alfven surface lying farther away from the planet than the 
sonic one. 

Accretion streams onto the stellar surface can form also in 
cases where the planetary outflow remains weak and the stellar 
wind is stopped by the planetary magnetic field. This is illus¬ 
trated in Fig. [TTJ which shows the final configuration of model 
FVRB. This case is similar to model FVRb, shown in the right 
panel of Fig. 0 in that a combination of high UV flux and large 
escape speed generate a dense discharge from the planetary sur¬ 
face but only a weak outflow. In this case, however, the stronger 
magnetic field causes the interface between the planetary mag¬ 
netosphere and the shocked stellar wind to lie at a larger distance 
from the planet’s surface; in particular, it now lies beyond the L\ 
point. At that location, the gravitational pull from the star and the 
thermal pressure of the dead zone combine to deform the plan¬ 
etary magnetic field lines in such a way that accretion streams 
are formed. These streams, however, do not flow smoothly: as is 
seen in the left panel of Fig. \TT\ their interaction with the stellar 
wind and the magnetic field of the star causes their trajectories 
to split several times before they reach the stellar surface. 

A similar situation characterizes our high-flux “near” mod¬ 
els, as illustrated in Fig. [12] with model FvrB. In this case, the 
planetary outflow is massive enough and the L\ point is located 
close enough to the planet’s surface that both a planetary tail 
and an accretion stream are formed, resembling the behavior of 
model FVRB. However, the different position of the planet rela- 
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Fig. 9. Same as Fig.[6l but for model FvRb. 
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Fig. 12. Same as Fig. [6) but for model FvrB. 


4. Classification of star-planet interactions 

Our simulations of the dynamical interaction between a mag¬ 
netized, wind-driving host star and a magnetized hot Jupiter 
that loses mass to photoevaporation were described in Sect. [3] 
in terms of the resulting morphological structures. We now at¬ 
tempt to distill from this phenomenology a general classification 
scheme that is based on the underlying dynamical processes op¬ 
erating in such systems. We quantify the relative influence of 
these processes with the help of characteristic length scales, and 
use the latter to identify four basic types of interaction. We then 
apply this scheme to the categorization of the models listed in 
Table [2 


five to the star modifies the evolution in this case. First, the frag¬ 
mented tail is not pushed away by the stellar outflow, which is 
still weak at this location. Instead, the stellar outflow mixes with 
these fragments and slows them down, leading to their accretion 
by the star. Second, the denser ambient stellar gas in this model 
also enhances the mixing with the accretion stream, slowing the 
latter down and causing it to hit the stellar surface close to the 
orbital location of the planet. Finally, as the orbital motion of 
the planet is now faster than the stellar rotation, the stream does 
not attain a steady-state configuration. Instead, the magnetic field 
topology undergoes a continuous readjustment, with new mag¬ 
netic accretion channels forming periodically. Model FVrB (not 
plotted) exhibits a similar morphology to that of model FvrB, 
with accretion flows onto the star originating from both the up¬ 
stream and downstream regions of the planet, although the larger 
escape speed results in less mass leaving the planet in this case. 
We however find that there is no transfer of mass to the star when 
the UV flux is low, irrespective of whether the escape speed is 
large or small (models fVrB and fvrB, respectively; these are 
also not plotted). 


4.1. Characteristic length scales 

The formation of a quasi-stationary morphological structure in 
the interaction between a star and a close-in planet entails the es¬ 
tablishment of pressure equilibrium between the stellar andplan- 
etary plasmas at the interface separating these two mediaQ The 
relevant frame of reference for considering this equilibrium is 
that of the planet, and we label the pressure on the stellar side of 
the boundary in this frame by F a mb- The characteristic distance 
of this surface from the center of the planet, measured along (or 
close to) the line to the stellar center, is determined by whether it 
is the magnetic pressure, outflow ram pressure, or thermal pres¬ 
sure that dominates on the planet’s side of the boundary. 

If the pressure of a closed planetary magnetosphere dom¬ 
inates, we can estimate the relevant scale (R m ) by assuming 


8 This interface formally constitutes either a contact or a tangential 
discontinuity. As confirmed by our simulations, such discontinuities 
may be subject to instabilities, which, among other effects, can induce 
mixing of the two plasmas. 
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a dipolar field, B 0 (R ) = B 0 (R 0 /R ) 3 , and equating its pressure 
Pi nago = B 2 JSn to / J a nib- This gives 


2 d6 \ 1 


R* 


( B 2 X \ 

\ &7tP amb / 


(28) 


If, however, the planetary outflow is strong enough to become 
supercritical by the time it is stopped (in a shock) by collid¬ 
ing with the stellar gas, the relevant scale (R w ) is given, in¬ 
stead, by equating P am b to the ram pressure P m mo = p 0 v 2 0 . 
Expressing the density in terms of the planetary mass loss rate 
Mo = 4nR^p wo v wo , we get 


R* 


I MqVwO \ 1/2 

\ 4ti P am b / 


(29) 


Which of these pressure components (magnetic or ram) domi¬ 
nates the interaction with the stellar plasma is reflected in the 
relative magnitudes of the associated characteristic radii. 

One could in principle also consider the contribution of the 
thermal pressure P = Ipk^T/m v of the planetary atmosphere, 
but in practice its role is negligible because of its expected 
rapid decline with radius. In particular, under the isothermal 
approximation, R 0 (R') oc exp(R 0 /R') (see Eq. [IT]), which rep¬ 
resents a much faster drop than that of the magnetic pressure 
(^mag 00 (Ro/R') 6 ). In fact, the magnetic pressure likely already 
dominates the thermal pressure at the (subsonic) base of the 
planetary outflow, where, using Jupiter’s mass, radius, and mag¬ 
netic field (B 0 - 4 G), and the base density (p 0 - 10~ 14 g cm -3 ) 
and temperature (T Q — 10 4 K) inferred from our high-UV-flux 
model (see Fig. [4]), we estimate 




P o 

R mago 



(30) 


Note, however, that when the stellar wind is stopped by a su¬ 
percritical planetary wind, the plasma on the planet’s side of the 
contact discontinuity is dominated by its thermal pressure (hav¬ 
ing passed through a shock). 

The ambient pressure can be approximated by the sum of 
four contributions: 


R amb — E w * + P mag R ram + R ram orb > (31) 

which represent, respectively, the stellar plasma’s thermal, mag¬ 
netic, and “intrinsic” ram pressure components, and the ram 
pressure associated with the relative motion between the planet 
and the ambient (stellar) gas. In this approximation, we take the 
direction of the stellar wind to be purely radial and the direction 
of the orbital motion to be purely azimuthal in the lab frame. 
If the stellar wind is supercritical, a bow shock will form at the 
interaction site. In the limit v or b v w *, the apex of this shock 
will lie along the line connecting the centers of the star and the 
planet, and the shock surface will be symmetric about this line. 
However, given that the orbital speeds of hot Jupiters are typi¬ 
cally >100 km s -1 and are thus highly supersonic with respect to 
the ambient gas, such a shock will form even when the stellar 
wind is still subcritical when it is stopped. In this case, however, 
the shock will be oriented at a finite angle to the aforementioned 
line (see Eq.[37]below), and its apex will be displaced away from 
this line (in the direction of the planet’s motion)0For typical pa- 

9 Note that the speed of the orbital motion for these systems is of the 
order of the terminal speed of the stellar wind, so that, when the interac¬ 
tion occurs in the supercritical region of the stellar wind, the resulting 
bow shock will be significantly inclined with respect to the instanta¬ 
neous “line of centers.” 


rameters, we expect that the last two terms in Eq. (ED constitute 
the dominant contributions to Pamb- 

One other characteristic length scale for this problem is the 
tidal (or Hill) radius, 

/ M \ 1/3 

* isi) (32) 


which gives the distance of the L\ Lagrange point from the 
planet’s center. This distance is relevant to the question of 
whether the morphology of the planetary gas that interacts with 
the stellar plasma will be shaped only by the relative motion 
between these two media (through the effect of the ram pres¬ 
sure terms in Eq. [31]) or whether the stellar gravitational field 
will also play a role. It can be expected that at least some of 
the planet’s plasma will be accreted onto the stellar surface if 
R t < max ( R m , P w ), but that essentially all of the planet’s gas 
will remain confined to the vicinity of its orbital radius if this 
inequality is not satisfied. 

Before concluding this subsection, we note that the compo¬ 
nents of Pamb can be expressed in terms of the stellar parameters 
just as was done above for the planet’s pressure components. In 
particular, using M* = 47rP 2 p w *v w *, B*(R) = B*(R*/R 0 r b) 3 , and 
the isothermal-gas assumption, we can express the terms appear¬ 
ing on the right-hand side of Eq. (1371) in the form 


mag * 


2 P W* 

m p 

Bj! R, 
S/r \R otb 


MtkftTy, 

2nm„R°- u \ 


— p w *V w * — 


ram orb 


-'w* 17 or b 


4 ^ rb v w * 


(33) 

(34) 

(35) 

(36) 


Thus, the values of the three characteristic length scales ( R m , P w , 
and R t ) can be estimated from a given set of planetary (. R 0 , M 0 , 
M-o ? Bo), orbital (Porb? ^orb)? and stellar ( R *, 44*, 44*, 

r w *, B *) parameters. 


4.2. Types of interaction 

We base our classification of the flow structures that arise from 
the hydromagnetic interaction between a close-in giant planet 
and its host star on the relative magnitudes of the three char¬ 
acteristic length scales considered in Sect. 14.11 R m , P w , and R t . 
The four basic types of interaction that we identify in this way 
are sketched in Fig. [13] and discussed below. We note, however, 
that any given system may not fall exclusively under a single cat¬ 
egory. This could be due to the planet having an eccentric orbit, 
to the stellar wind being nonaxisymmetric, or to the presence of 
multipolar components in the stellar magnetic field, which could 
each lead to variations in the values of R m and R w on a timescale 
of -days. A stellar magnetic cycle would be likely to induce 
variations on timescales of -years in the properties of the stellar 
magnetic field, the stellar wind, and the stellar UV flux. Changes 
on even longer timescales could result from the effects of stellar 
evolution. 

Type I: R t > R m > R w 

Type I interactions occur when the planetary outflow is weak 
(corresponding to either a low irradiating flux or a large escape 
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Type I: bowshock and thin tail 


Type II: colliding winds and tail 




Type III: strong planetary-wind, 
accretion and tail 


Type 



IV: Roche-lobe overflow 
accretion and tail 


/ 



Fig. 13. Schematic of the different possible star-planet interactions, showing face-on views in the orbital plane of the four distinct 
morphological structures (denoted by the Latin numerals I, II, III, and IV) as they appear in the planet’s frame. The large and small 
solid disks (on the left and right sides of each panel) represent the star and the Hot Jupiter, respectively, the shaded areas highlight 
material that flowed out of the planet, the arrows indicate gas motions, and the closed loops sketch the planetary magnetosphere. 
The solid circular arcs in panels I and IV have a radius R m , the distance where the magneto spheric pressure equals the total ambient 
pressure, whereas the dashed circular arcs in panels II and III mark the distance R w where the ambient pressure equals the ram 
pressure of the planetary wind. The dotted line indicates the contour of the effective (gravity plus centrifugal) potential that crosses 
the Lagrange point L\ (at a distance R t from the planet). The proposed classification scheme is based on the relative ordering of R m , 
R w , and R t . 


speed), so that the stellar plasma is intercepted by the planetary 
magnetic field (i.e., R m > R w ). As we noted in Sect. 14.11 the rel¬ 
ative motion between the planet and the ambient gas is typically 
high enough to lead to the formation of a bow shock upstream of 
the planet. The shocked stellar flow sweeps back the planetary 
magnetic field and diverts the planetary outflow in the down¬ 
stream direction, leading to the formation of a thin planetary tail. 
The stronger the planetary field, the larger its pressure — hence 
the farther away from the planetary surface does the shock ap¬ 
pear ( R m is larger), and the wider the tail that is produced. The 
tail gas remains confined to the vicinity of the planet’s orbit be¬ 
cause R t > R m . 

The orientation of the tail is determined by the direction of 
the incident stellar wind as seen in the corotating frame at the 
apex of the bow shock. Approximating the wind velocity to be 
nearly radial in the lab frame and of magnitude v w (see Fig. 0), 
the angle that the tail makes with respect to the tangent to the 


planet’s orbit (indicated on the top left drawing in Fig. |T3]) can 
be approximated by 

n 4- V w* / V w*^orb'\ 

6 t aii ~ arctan-= arctan ——— . (37) 

V 0 rb \ GM * J 

(see also IVidotto etal.! l 2010h . In the limit of an extended and 
static stellar corona (v w * = 0), the tail will trail the trajectory, 
forming a torus around the star. In the other limit, when a strong 
wind hits a comparatively distant planet, the tail will be almost 
perpendicular to the orbit, and so will lie nearly along the line of 
sight during transits. Note that we did not include the effect of ra¬ 
diation pressure, which is the dominant force in the formati on of 
comet-type tails in close-in rocky planets (e.g., Rapp aport et al.l 
l2012ll2014l) . In the case of the gaseous giant planets considered 
in this paper, radiation pressure on Lya lines could potentially 
play a role, but its contribution would be negligible if (as we 
suggest below) the gas in the tail is highly ionized. Furthermore, 


13 


















T. Matsakos et al.: Star-planet interactions 


even in regions where the gas may still be only partially ion¬ 
ized (for example, near the base of the planetary outflow), the 
amount of gas affected b v this process would be stro ngly limited 
by opacity effects (e.g.. iTremblin & Chian gl [2013b . It is, how¬ 
ever, conceivable that radiation pressure on resonance lines such 
as those of H I and Mg I could contribute to the blueshifted 
spectral features observed in the meas ured a bsor ption pro- 
f iles of these lines (e.g., iBourrier & Lecavelier des Etangsll2013l: 
iBourrier et al .112014 and references therein). 

Our simulations indicate that the gas that dominates the col¬ 
umn density between the Hot Jupiter and the bow shock is of 
planetary origin, rather than the shocked stellar plasma. We ob¬ 
tain a representative value for the hydrogen column density of 
this component of Aho ~ 8 x 10 15 cm -2 (using a typical hydro¬ 
gen number density of «ho ~ 6 x 10 5 cm -3 and a path length 
of ~ 2R Q , based on the simulation results for model fVRb; see 
Fig. [8]). This column can be compared with the equilibrium col¬ 
umn of an ionized slab with this density, 


Nmi 


O'BftHoiqjV 


(38) 


9x 10 19 


/ F l ,v 

\( V 

\ 450 erg cm -2 s -1 / 

/ \ 6 x 10 5 cm -3 / 


cm 2 , 


where - 2.6 x 10“ 13 cm 3 s -1 is the Case B recombina tion 
coefficient, and where (following iMurrav-Clav et al1l2009h we 
adopted a characteristic ionizing photon energy of £uv = 20 eV. 
This estimate indicates that, even for a comparatively low ioniz¬ 
ing flux, the entire planetary gas in the interaction region would 
be fully ionized, implying that this component would not be de¬ 
tectable in Ly a. A possible caveat to this conclusion may arise 
when the configuration is not in quasi-static equilibrium, and in 
particular if the ionization time, 


fionize 


/ £ uv \ / £uv \ 3 
\cr 0 Fuv/ \13.6eV7 

3 g - 10 4 f - — 

\20eV/ \ 450 erg cm -2 s _1 


s 


(39) 


(where <xo - 6x 10 -18 cm 2 is the hydrogen photoionization cross 
section at the Lyman edge), is longer than the travel time of the 
planetary gas to the interaction region. However, this is probably 
not relevant to Type I configurations on account of the low ve¬ 
locity of the outflowing gas. The planetary gas in Type I systems 
might nevertheless be detectable in other absorption lines, such 
as the UV h and k resonance lines of Mg II. 

Type II: R t > R w > R m 


Type II interactions occur when the planetary outflow is com¬ 
paratively strong (corresponding to either a high irradiating flux 
or a small escape speed), so that the stellar plasma is intercepted 
by the planetary outflow (i.e., R w > R m ). We again expect this 
plasma (which moves with a velocity of a few hundred km s -1 
with respect to the planet) to pass through a bow shock, and, after 
being shocked, to sweep back the planetary outflow into a tail. 
Because of the larger momentum flux in the planetary outflow 
in comparison with the Type I case, Type II tails are wider than 
their Type I counterparts for a given value of B 0 . However, just 
as in the case of Type I interactions, the tail gas remains confined 
to the vicinity of the planet’s orbit, in this case because R t > R w . 
Its tilt angle (ftaii) with respect to the orbit can be similarly esti¬ 
mated using Eq. (1771) . 

Since in this case the planetary outflow is supercritical by 
the time it collides with the stellar plasma on the day side of the 


planet (at which point its speed v wo is < 30kms _1 , on the or¬ 
der of v esco ), it is also decelerated in a shock. The two shocks 
are separated by a contact discontinuity, across which the total 
(thermal + magnetic) pressure is the same. This, in turn, implies 
that the normal components (with respect to the contact discon¬ 
tinuity surface) of the ram pressures of the two flows are approx¬ 
imately equal in the frame of the contact discontinuity (which, 
in a steady state, is the same as that of the planet). As a rough 
estimate, we write 

Pw*(v w * + v or k) ~ p W oV wo , (40) 

where, as we already noted, v w * is typically of the order of v or b 
for a Hot Jupiter. The estimate (l40l) implies that the ratio of 
the preshock densities in the interaction region is p w */p W o ~ 
v* 0 /(Vw, + v 2 rb ), which for v wo ~ 0.1 v w , is -0.01. The ratio 
of the corresponding postshock densities need not be the same 
on account of differences in the shock strength and the post¬ 
shock cooling efficiency between the two winds, but it would 
likely still be 1. It also follows from Eq. (l40l) that the ratio 
of the mass inflow rates into the interaction zone, ~ p w *(v 2 * + 
v^r b ) 1/2 /(PwoVwo), is - v wo /04* + v orb) 1/2 V" 0 - 1 for our adopted 
values). This suggests that, as in the Type I case, the postshock 
density and column density in the interaction region is domi¬ 
nated by the planetary material. Although our simulations have 
not produced examples of a Type II interaction, we can use the 
simulation results for model FvRb (Fig.0, which represents an 
example of a Type III interaction that involves a comparatively 
massive planetary outflow (see Table [2] and Sect. 14.3b , to infer 
an upper bound of -10“ 16 g cm -3 on the typical postshock plan¬ 
etary density in this case. We note in this connection that a su¬ 
percritical planetary outflow is also shocked on the night side of 
the planet when it hits the gas that had been shocked further up¬ 
stream and subsequently dragged backward around the planet0 
Moving farther out from the star in that direction, one again en¬ 
counters unshocked stellar-wind gas that expands nearly radially, 
so the density there rapidly becomes very low. 

Our estimate of the density of the planetary gas in the 
interaction region is two orders of magnitude larger in this 
case than the typical value we adopted for the Type I inter¬ 
action (~10 -16 gem -3 , associated with a high flux value, vs. 
~10 -18 gcm -3 , which was obtained for a low UV flux). Using 
Eq. (l38b , we deduce a nominal ionized column of Nmi ~ 
10 21 cm -2 for these parameters, as compared with a total col¬ 
umn of Aho < 8 x 10 17 , indicating that, in this case too, the 
interacting gas can be expected to be fully ionized if in equilib¬ 
rium. The equilibrium assumption could be questioned since the 
gas is clearly in motion for Type II configurations; however, for 
the high flux value that characterizes the flow in model FvRb, 
the nominal ionization time inferred from Eq. (l39b (a few tens of 
seconds) is much shorter than any relevant flow time. In addition, 
the shock heating of the planetary outflow can also contribute to 
the ionization of this gasQ It is therefore likely that much of the 
gas that participates in a Type II interaction is ionized, although 
this conclusion needs to be checked explicitly for each given set 
of parameters. 


10 The detailed morphology of the night-side region may change if 
the assumption of equivalent outflows from the two sides of the planet 
is modified. However, we do not expect our basic classification scheme 
to be sensitive to these details. 

11 The postshock temperature scales as the square of the upstream gas 
velocity normal to the shock front; one can therefore expect tempera¬ 
tures of a few times 10 4 K for the shocked planetary flow, as compared 
with a few million degrees for the shocked stellar plasma. 
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Type III: R w > R m & R w > R t 

Type III interactions are similar to those of Type II in that they 
involve a strong planetary wind ( R w > R m ). However, in this 
case the planetary outflow, after being stopped in a shock, is only 
partly swept back into a tail, with another part responding to the 
gravitational pull of the star (which is felt because R w > R t ) 
and penetrating the stellar plasma ahead of the planet, forming a 
stream that spirals in and accretes onto the star0 

The accretion stream exhibits a complex, fragmented mor¬ 
phology that arises from Kelvin-Helmholtz and Rayleigh-Taylor 
instabilities — triggered, respectively, by the tangential veloc¬ 
ity shear and by the radial acceleration exerted by the stellar 
wind gas on the denser planetary material. The stellar gravity 
further contributes to the development of the latter instability 
on account of the “heavy on top of light” density stratification. 
These processes, and possibly also the magnetic stresses acting 
on the inflowing plasma, fragment the inspiraling stream, result¬ 
ing in the formation of multiple accretion filaments that hit the 
stellar surface at different spots. The general location where the 
accretion stream impacts the stellar surface is determined by the 
stellar and planetary wind parameters and by the orbital distance 
of the planet, but it is typically well ahead of the instantaneous 
sub-planetary point (the phase difference being -90° in the ex¬ 
ample shown in Fig. 0. We do not observe the formation of a 
circumstellar disk in our simulations, evidently because the drag 
exerted by the stellar outflow and the mixing between the two 
plasmas result in efficient removal of angular momentum from 
the inspiraling gas. 

There have been several claims in the literature for enhanced 
activity on the host star — manifested in optical chromospheric 
emission lines (particularly Ca II), transition-region FUV emis¬ 
sion lines (such as Si IV), and coronal X-ray emission — that 
migh t be attributable to an interaction with a close-in planet 


(e.g., Shkolni 

ket al. 2003,2004, 

200l 20Q8j;[Walker et al. 2008; 

Pillitteri et al. 

201(1120111 2014 

, 2015k The accretion streams 


produced in Type III interactions could potentially trigger such 
events. In particular, they may be associated with active regions 
that are inferred to lie significantly ahead of the sub-planetary 
point — as in HD 179949 and r Boo, where a phase difference 
of -70° was inferred from measurement s of the Ca II K line 
dShkolnik et al .1120081: IWalker et al.1120081) . and in HD 189733, 
where a lead phas e of ~70°-90° was deduce d from X-ray and 
FUV observations dPillitteri et al .12014120151) . An alternative in¬ 
terpretation of these phenomena that has been discussed in the 
literature involves a direct magnetic interaction of the type ob¬ 
served in planetary moons (such as Io) in the solar system. Note, 
however, that in their simplest form, such interactions are pre¬ 
dicted to occur near the sub-planetary point. Furthermore, the 
energy flux generated in such an interaction in a source like 
HD 179499 is calculated dSaur et al.ll2013l) to be over two or- 
d ers of magnitude lower than the observationally inferred value 
dShkolnik etal.ll2005h . 

Type IV: R m > R w & R m > R t 

This case is analogous to the Type III one in that the interaction 
region lies outside the planet’s Hill radius, resulting in the ac¬ 
cretion of planetary gas onto the stellar surface. However, just 

12 If the planet orbits close enough to the critical surface of the stellar 
wind (as defined in the corotating frame), it may be possible for the 
stellar outflow to encounter the boundary of the stream while it is still 
subcritical, in which case no bow shock will form. However, a shock 
may still form when the stellar wind intercepts the planetary tail if it is 
already supercritical at that location (see, e.g., Fig. 0. 


as a Type III interaction could be regarded as a special case of 
a Type II one, with both involving a strong planetary outflow, a 
Type IV interaction also bears similarity to a Type-I interaction 
in having a comparatively weak planetary wind ( R m > R w ). In 
this case the planetary magnetosphere, which is loaded by the 
subsonic gas evaporated by the stellar radiation field, intercepts 
the stellar plasma beyond the L\ point, where the stellar grav¬ 
ity is the dominant force acting on the planetary gas. This situa¬ 
tion resembles the classic Roche-lobe overflow picture, although 
both the planetary magnetic field and the outflowing, magnetized 
stellar plasma also play a role in shaping the accretion flow mor¬ 
phology. In particular, when the orbital radius is small enough, 
the relative motion between the stellar gas and the planet may be 
supercritical and therefore lead to the formation of a bow shock 
ahead of the planet, as in the Type I case, even as mass transfer 
from the planet to the star takes place in the vicinity of the sub- 
stellar region. However, no bow shock forms on the star-facing 
side in our “near” models because the stellar magnetic field (and 
hence the Alfven speed) at R m is large enough in these cases to 
keep the relative motion between the stellar gas and the planet 
subcritical (see Fig. m. 

The structure of the accretion flow and its dynamical impli¬ 
cations are similar in this case to the situation in Type III inter¬ 
actions. However, because of the absence of a strong planetary 
outflow in the Type IV case, the accretion stream does not exhibit 
the spiral structure observed in Type III interactions and instead 
maintains a nearly linear shape, hitting the stellar surface near 
the sub-planetary point. This type of interaction may thus be rel¬ 
evant to the interpretation of enhanced stellar chromospheric/X- 
ray activity that is inferred to occur near the sub-planetary point 
(and could provide an alternative to the “direct magnetic inter¬ 
action” scenario also in this case). 

A Roche-lobe overflow model was previously proposed for 
the interpretation of early ingress indications in UV absorption¬ 
line observations o f the short-period (R or b = 0.023 AU) Hot 
Jupiter WASP-12b dFossati et al.ll2010h . In particular. lLai et al.1 
(120101) . considered the possible contributions of an accretion 
stream and of an inner accretion disk, and also speculated on 
absorption in the interaction region between a stellar wind and 
the planetary magnetosphere. They did not, however, include the 
effect of the stellar radiation heating of the planet and the influ¬ 
ence of the stellar plasma on the accretion stream. iBisikalo et aD 
( 2013 ) carried out a hydrodynamical simulation and drew atten¬ 
tion to the interaction between the accretion stream and the stel¬ 
lar outflow, but they did not consider the effect of a stellar and/or 
a planetary magnetic field. The results presented in this paper 
provide a more general framework for modeling systems like 
WASP-12, and may also aid in modeling other types of close- 
in planets where Roche-lobe overflow could play a role (e.g., 
IValsecchi et al.lll2014l) . 


4.3. Classification of the simulations 

The application of our proposed classification scheme to the 
models that we have simulated is presented in Table [3] Note 
that Type II interactions are not represented in this table. This 
interaction is intermediate between the Type I and Type III cases 
— for which we do have examples — in that the planetary out¬ 
flow must be strong enough to cause R w to exceed R m but not 
so strong that it will also exceed R t . Computational constraints 
prevented us from exploring a wider parameter space that would 
have encompassed this interaction. 
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Table 3. Classification of the models listed in Table [2] 


Model 

Interaction type 

FvRB 

III 

FvRb 

III 

FvrB 

IV 

fvRB 

I 

fvRb 

I 

fvrB 

I 

FVRB 

IV 

FVRb 

I 

FVrB 

IV 

fVRB 

I 

fVRb 

I 

fVrB 

I 


5. Summary & conclusions 

In this paper we study the different types of star-planet interac¬ 
tions that can occur in systems harboring Hot Jupiters. We per¬ 
form a series of parametrized 3D MHD numerical simulations, 
incorporating the star, the planet, and consistent stellar and plan¬ 
etary winds. Based on the results of a grid of models, we propose 
a classification of the interactions into four general categories 
(which we label I-IV; see Fig. |T3]). We describe in detail the dy¬ 
namical features exhibited in each case, and provide expressions 
for the three characteristic length scales that underlie our classi¬ 
fication scheme in terms of the system’s physical parameters. We 
also briefly discuss general characteristics of these interactions 
that could have a bearing on their observational properties. 

The main conclusions of our analysis can be summarized as 
follows: 

1. Except very close to the star, where the stellar magnetic field 
(and the associated Alfven speed) might be large enough to 
render the flow subcritical, the stellar outflow intercepted by 
the planet is stopped in a bow shock. The shock forms even 
when the stellar wind is not yet fully accelerated because of 
the high (> 100 km s _1 ) azimuthal velocity component of the 
stellar gas in the planet’s frame. 

2. In Type I interactions, the stellar plasma is intercepted by the 
planetary magnetosphere, which, in turn, is swept back into 
a plasma tail that is blown away by the stellar wind. 

3. In interactions of Type II and III, the stellar plasma is in¬ 
tercepted by a strong planetary outflow that opens up the 
magnetosphere. The shocked planetary plasma is dragged 
backward by the stellar outflow, forming a wide tail that is 
fragmented by dynamical instabilities. In the Type II case, 
the shocked planetary plasma lies inside the Roche lobe and 
remains confined to the vicinity of the orbital radius. In con¬ 
trast, in the Type III case the interaction occurs outside the 
Roche surface, resulting in part of the shocked planetary 
plasma forming an accretion stream that spirals toward the 
star until it impacts the stellar surface. The dynamical inter¬ 
action of this stream with the magnetized stellar wind en¬ 
ables this gas to lose energy and angular momentum and fall 
in, instead of forming a plasma torus. The stream is also sub¬ 
ject to dynamical instabilities, which cause it to split into 
multiple accretion filaments that hit the stellar surface at dif¬ 
ferent spots. 

4. In Type IV interactions the stellar outflow is intercepted by 
the planetary magnetic field outside the Roche radius. This 
situation resembles a classic Roche-lobe overflow, although 
it differs in detail on account of the presence of the planetary 
magnetic field and of the magnetized stellar outflow. In con¬ 


trast with the Type III accretion stream, which reaches the 
stellar surface well ahead of the sub-planetary point, in this 
case the planetary gas falls in nearly radially. 

5. In general, the density and column density in the interaction 
region are dominated by the planetary gas. Typically, this gas 
can be expected to be fully ionized, implying that detection 
through absorption features of neutral species such as H I 
and Mg I is unlikely. 
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